Observations on Sound Propagation in Rapidly Rotating Bose-Einstein Condensates 
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Repulsive laser potential pulses applied to vortex lattices of rapidly rotating Bose-Einstein con- 
densates create propagating density waves which we have observed experimentally and modeled 
computationally to high accuracy. We have observed a rich variety of dynamical phenomena rang- 
ing from interference effects and shock-wave formation to anisotropic sound propagation. 
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Introduction. — Supcrfluids can be described in terms of 
an order parameter associated with the superfluid phase 
of the system. In the case of dilute Bose-Einstein con- 
densates a natural choice for the order parameter is the 
macroscopic wavefunction, whose squared modulus yields 
the particle density of the condensate. It follows that the 
superflow, given by the gradient of the phase of the con- 
densate, must be irrotational and therefore rotation of 
a superfluid requires the presence of quantized vortices 
with cores of vanishing density Q. The minimum en- 
ergy configuration of a rotating superfluid is a triangular 
Abrikosov lattice of vortices as was originally established 
in the context of Hc-II. This fundamental superfluid be- 
havior has also been witnessed in dilute Bose-Einstein 
condensates 0HSE1. 

The structure and dynamics of single vortices in dilute 
Bose-Einstein condensates have been extensively studied 
0. The presence of multiple vortices radically changes 
the spectral properties of the condensate while a rela- 
tive motion of the interacting vortices enriches the su- 
perfluid dynamics. Recently collective behavior of vor- 
tex lattices has attracted much theoretical and exper- 
imental interest. Oscillations of the giant vortex core 
of vibrational Tkachenko modes 
and the structure of the con- 
densate in the lowest Landau level 0, 0, 0] have been 
observed by changing the moment of inertia of the rotat- 
ing condensate by a selective atom removal. The tilting 
mode of a vortex array has been observed by applying a 
method analogous to the driving of the scissors mode of 
a non-rotating condensate [lgj . 

In dilute Bose-Einstein condensates, compressional 
sound waves correspond to modulations in the conden- 
sate density. Propagation of such sound waves were ob- 
served in harmonically trapped ground state condensates 
[lf| . and in two-component condensates density defects 
were seen to deform into vortex rings and dark solitons 
[ifl l2l1 | . The effect of an optical lattice on sound prop- 
agation has also attracted attention recently Since 
the condensate density vanishes at vortex cores, the rota- 
tion of the condensate is expected to have notable effects 
on the sound propagation. In this Letter we study both 




experimentally and theoretically the response of a rapidly 
rotating Bose-Einstein condensate to a variety of exter- 
nal repulsive laser pulses. We compare our experimental 
observations and a priori time-dependent numerical sim- 
ulations, and show how the topological constraints set by 
the vortices result into exotic non-equilibrium dynamics 
such as anisotropically propagating sound waves. 

Experiment and theory. — The experimental results 
presented in this paper are obtained by starting with 
rapidly rotating 87 Rb condensates held in a cylindrically 
symmetric trap with trapping frequencies {ujj_ , u> z } = 
2ir x {8.3, 5.3} Hz. The condensates typically contain 
N = 3.5 x 10 6 atoms and are rotating around the z-axis 
with a rotation rate of 0.95 uj±. To create blastwaves, 
a far off-resonant laser beam, centered on the conden- 
sate, is sent through it along the z-axis. This beam has 
a wavelength of 660 nm, which is far detuned from the 
Rb D2 line at 780 nm. It creates a dipole potential given 
by 
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such that for a beam with power P = 0.64 mW and waist 
d = 19 (im as used for Fig. the atoms at the center of 
the laser beam experience a repulsive potential with char- 
acteristic frequency of 2tt x 78 Hz. The potential height 
in this case surpasses the chemical potential of the con- 
densate by a factor 39. After the pulse, the condensate 
is expanded by an anti-trapping geometry as described 
in |23j, and imaged along the axis of rotation. 

The experiment is simulated in the mean-field formal- 
ism by calculating the full 3D evolution of the condensate 
wavefunction with the time-dependent Gross-Pitaevskii 
equation. We adopt fko± for the energy scale and oq = 
yjh/2muji_ for the spatial length scale. From the compu- 
tational point of view, in the absence of the laser field, 
the condensate is fully characterized by the dimension- 
less values for the nonlinearity constant C = gN/huj±aQ, 
and the trap anisotropy parameter A = lo z /uj^, when the 
wavefunction ip is normalized to unity. Above, g is the 
usual s-wave coupling constant. Modeling this experi- 
ment we choose C = 1.6 x 10 5 and A = 0.64. 



Wide pulse. — Figure ^ shows experimental observa- 
tions (upper row) and their theoretical counterparts 
(lower row) of the integrated condensate density, for the 
case where a laser pulse of power P = 4.5 mW and waist 
d = 60 //in (roughly the radius of the condensate) is ap- 
plied for varying lengths of time. The pulse durations 
are 5 ms (a,d), 10 ms (b,e), and 30 ms (c,f). For the 
shortest pulse (a,d), mainly the collective modes such as 
the Tkachenko waves and longitudinal breathing modes 
are excited. When the duration of the pulse is doubled 
(b,e), the central part of the condensate is repelled to 
form a narrow outer band of high density where the con- 
densate is channeled between the vortex planes forming 
supcrfluid flares. The vortices are compressed both in 
core size and spacing there. In contrast, vortices in the 
central region experience an outward movement and ex- 
pand in both core size and spacing due to the strongly re- 
duced condensate density. It is notable that the number 
of particles per vortex in the centermost region remain 
practically unaltered from its equilibrium value. For even 
longer pulse duration (c,f) a five-fold increase in the trap- 
ping energy of the condensate takes place. When this ex- 
cess energy is converted into kinetics of the condensate, 
strong radial currents result and the stripes seen in (c,f) 
emerge due to the fact that the vortex cores constitute 
forbidden positions for the supcrfluid flow. The outer 
higher density halo at the condensate horizon is a rem- 
nant of the earlier state when the condensate assumed a 
shape of a thin cylinder. 

The physics seen in Fig. [2 is conveniently analysed us- 
ing Fig. [3 where different components of energy are plot- 
ted as functions of time corresponding to the simulation 
in Fig. IHd). All of the three cases, Fig. [TJd)-(f), ex- 
hibit similar charactristics in terms of absorbtion of en- 
ergy from the beam followed by a subsequent oscillatory 
exchange of kinetic and potential energies at the char- 
acteristic radial compressional mode frequency of 2.0uj± 
|24j . However, the amount of energy absorbed by the con- 
densate from the beam depends on the pulse duration, 
resulting in the qualitatively different features between 
the tree different cases in Fig. ^ 

The anti-trapping stage is not modeled in the simu- 
lations presented due to the limitation of computational 
resources. This may slightly alter the apparent size of 
the vortices due to the axial evolution as explained in 
[23| . However, as is shown in Fig. [21 the mean-field 
energy amounts to only a few percent of the total en- 
ergy and thus the dynamical evolution is mainly gov- 
erned by the interplay between kinetic and potential en- 
ergies, which explains the similarity between the in-trap- 
and anti-trapped results, apart from the scaling in space 
and time. 

Narrow pulse. — When applying laser pulses with nar- 
row beam waist (with respect to the width of the con- 
densate), travelling density perturbations are created. 
Strong applied power produces shock- waves, fringes and 
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FIG. 1: Response of a vortex lattice to wide repulsive laser 
pulses. Upper row displays experimental pictures and the cor- 
responding frames below them are from 3D numerical simu- 
lations. The pulse durations are 5 ms (a,d), 10 ms (b,e), 
and 30 ms (c,f). The experimental (computational) pictures 
are taken after 45 ms of anti-trapped expansion (respectively, 
(12, 16, and 22) ms in-trap evolution) after the end of the 
pulse. 
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FIG. 2: Kinetic Eu, trap Eti, and mean-field E m { energies and 
their sum as functions of time. The data is obtained from the 
same simulation as Fig.0^d). 

stripes as highlighted in Fig. [21 whereas if a weak pulse 
is used, smoothly propagating sound waves are formed 
as demonstrated at the end of the paper. The physics in 
Fig. is again characterised by the energy development. 
At early times, the kinetic energy is greatly increased 
as the strong laser potential pushes atoms at supersonic 
speeds away from the beam focus thereby generating a 
steep wavefront that rapidly propagates radially while 
experiencing a twist under the influence of the Corio- 
lis effect. The formation of the shock front at the edge 
of the beam is manifested by a sudden increase in the 
mean-field energy shortly after the end of the pulse as 
shown in Fig. |2{c) . The average initial velocity (v) of the 
repelled atoms may be crudely estimated by distributing 
the kinetic energy absorbed from the beam AEk over the 
particles within the beam volume V = 27rd 2 zxF- The re- 
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FIG. 3: Experimental (a) and theoretical (b,c) response of a 
vortex lattice to a localized 5 ms long repulsive laser pulse. 
Beam parameters are P = 0.64 mW and d = 19 /im. The 
experimental image is taken after 45 ms anti-trapped expan- 
sion, whereas the 3D computational image is taken after 14 
ms of in-trap evolution. 



suit is (v) = c s y / 2Cfiuj± AE^a^/V/i 2 , where c s = y/ \xjm 
is the speed of sound. This provides a prediction for the 
shock formation and for the case in Fig. ^a) (v) w 0.7c s 
whereas for Fig.[3]we obtain (v) ~ 3.7c s . After the shock 
is formed in Fig. |3J the hole created in the center of the 
condensate starts to fill in due to the refocusing effect of 
the trap. Meanwhile, the original outward propagating 
wavefront has broken into an outer layer of supersonic 
ripples and an inner primary wavefront moving at the lo- 
cal speed of sound, see also Fig. In the computational 
frames of Figs. Q] and the laser potential is off-centered 
3 /zm from the trap/lattice symmetry axis. 

Sound propagation. — In order to isolate the effect of 
vortices on the propagation of density modulations, we 
simulate perturbations in both the ground state conden- 
sate and the rotating vortex lattice and compare the re- 
sults. In Fig. 0] large amplitude density waves are gener- 
ated using a 5 ms long pulse with 0.15 mW power and 
13 /im beam waist. The subsequent time-evolution of 
the condensate density in the two cases can be compared 
in the sequence of frames for the ground state (lhs) and 
the vortex lattice state (rhs). The frames for the ground 
state are taken along the line y = z = 0, whereas for 
the lattice the densities in the plane z = are averaged 
over circles of radii r in order to integrate out the large 
density fluctuations due to the vortex cores. To facilitate 
comparison, the ground state and the vortex lattice state 
have the same number of atoms. However, the radial trap 



frequency for the ground state has been weakened by the 
factor y/l — (fi/wjj 2 in order to compensate for the cen- 
trifugal effect experienced by the vortex lattice. 

In both cases a shock front with a width of the order 
of the healing length of the condensate is generated, as 
displayed in the second frames at t = 4 ms. As this wave 
reaches the boundary of the driving potential it suddenly 
slows down to the local speed of sound emitting a col- 
lection of high frequency excitations travelling in front 
of the primary wavefront. These supersonic ripples carry 
away part of the kinetic energy loaded in the shock-wave. 
Similar density modulations at the front edge have been 
recently reported [3,|2i| and were speculated to be a fail- 
ure of the mean-field theory j^j. However, these fringes 
are clearly present also in the experiments, as is shown in 
Fig. |3J and are a real physical effect. After a quarter of 
a trap period, the nonrotating condensate forms a den- 
sity peak in the vicinity of the trap center accompanied 
by sharp fringes. In the vortex lattice case, however, a 
centrifugal barrier exists and therefore the center of the 
trap remains nearly particle free at this stage of the evo- 
lution. Similar fringes to those seen in the ground state 
are also present in the lattice, but they are not visible 
in the (rhs) of Fig. 01 because of the above mentioned 
averaging. They can, however, be seen in Fig. 

In the vortex lattice, subjected to the same small am- 
plitude perturbations as the ground state, we find that 
the average radial speed of sound is equally well described 
by the equation for c s evaluated using the chemical po- 
tential of the lattice state in place of the one for the 
ground state. Alternatively, the influence of vortex lat- 
tice on the speed of sound may be understood via renor- 
malization of the effective interaction energy by a factor 
of b > 1. The value of b (and the chemical potential) de- 
pends on the angular velocity of the condensate, resulting 
in b w 1.16 in the quantum Hall regime ^3] an d hence to 
a faster speed of sound with respect to a ground state of 
the same spatial extent and particle number. This pre- 
diction is also in accordance with the velocities measured 
from our simulations and with the fact that the chemical 
potential of the lattice state is a factor of 1.1 larger than 
that for the corresponding ground state. 

Anisotropic sound. — An interesting feature of sound 
wave propagation in the vortex lattice is the appearance 
of an anisotropic wavefront. Initially, the sound wave re- 
flects the circular symmetry of the laser beam, but each 
time the wavefront encounters an orbital of vortices, it 
deforms to a hexagonal shape following the lattice planes, 
as depicted in Fig. where both the experimental and 
computational results are shown. On passing through the 
hexagonal grid of vortices, the wavefront becomes more 
circular, but is then deformed back to hexagonal on en- 
countering the next ring of vortices. We note that the 
supersonic ripples propagate more isotropically, indicat- 
ing that the lattice structure is relatively transparent to 
these high frequency modulations. 



4 





ms 






4 ms 






8 ms 






13 ms 


>\ /V 




17 ms 






21 ms 






25 ms 






30 ms 






34 ms 





-80 -40 40 80 -80 -40 40 80 

r [/im] r \fim] 



FIG. 4: Simulations (3D) on sound propagation in a ground 
state (lhs) and in a vortex lattice (rhs). The ordinates are 
normalized according to the temporal peak density. 



ical simulations. Although mean-field theory is known 
to describe well dilute Bose-Einstein condensates at low 
temperatures, it is perhaps surprising how accurately 
it predicts the behavior of rapidly rotating condensates 
even in the presence of shock waves and after the lattice 
structure seems to have been completely destroyed. We 
have identified a sudden increase in the mean-field en- 
ergy as a signature of the formation of a shock front and 
obtained a quantitative criterion for it in terms of the 
kinetic energy absorbed from the beam. The propaga- 
tion velocity of the peak condensate density is found to 
be well described by the Bogoliubov speed of sound for 
both the lattice and the ground state even in the case of 
large amplitude perturbations. Finally, we have discov- 
ered anisotropically propagating sound waves in vortex 
lattices and showed their inherent relation to the super- 
flow. 

Financial support from the Academy of Finland, and 
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FIG. 5: Anisotropic sound propagation as observed in an ex- 
periment (a) or in the 3D simulation (b). The duration of the 
pulse with P — 300 (jW, d = 18/xm is 2.5 ms. Experimen- 
tal (computational) image taken after 60 ms of anti-trapped 
expansion (16 ms in-trap evolution). 



The observed anisotropic sound propagation cannot be 
explained purely in terms of the static spatial density 
fluctuations due to the vortex cores. We have shown this 
by carrying out a simulation for a ground state conden- 
sate with a rigid grid of pinning potentials which mimic 
the density effect of a vortex lattice without phase circu- 
lation. In this case the sound waves propagate isotropi- 
cally without any sign of hexagonal deformations. Thus 
the development of hexagonal wavefront depends in an 
essential way on the quantum mechanical current around 
the vortex cores, the details of which will be discussed 
elsewhere. 

Discussion. — We have experimentally studied non- 
equilibrium dynamics of rapidly rotating Bose-Einstein 
condensates, triggered by repulsive laser pulses. These 
experiments arc modeled to high accuracy using numer- 
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